esaligna=read.table("esaligna.csv", head=T, sep=",")

par(mfrow = c(1,2))
boxplot(esaligna$dap~esaligna$talhao, xlab = "Talho", ylab = "DAP da rvore")
text(x = 1, y = 22 , "a")
dap_talhao <- matrix(data = c(esaligna$dap, esaligna$talhao), ncol = 2, 
                     nrow = length(esaligna$dap), byrow = FALSE)
med.dap <- c(tapply(esaligna$dap, INDEX = list(esaligna$talhao), FUN = mean))
sd.dap <- c(tapply(esaligna$dap, INDEX = list(esaligna$talhao), FUN = sd))

plot(x = NULL, y = NULL,  xlim = c(0.5,6.5), ylim = c(min(med.dap)-4, max(med.dap)+8), 
     type = "n", xaxt = "n", ylab = "DAP da rvore", xlab = "Talho")
axis(side = 1, at = c(1:6), labels = c("16", "17", "18", "22", "23", "32"))
text(x = 1.3 , y = 22.1 , "b")
points(med.dap)
segments(x0 = (1:6), x1 = (1:6), y0 = (med.dap[1:6]-sd.dap[1:6]),
         y1 = (med.dap[1:6]+sd.dap[1:6]))
